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ABSTRACT 

In this work we investigate the problem concerning the presence of additional bod- 
ies gravitationally bounded with the WASP-3 system. We present eight new transits 
of this planet gathered between May 2009 and September 2011 by using the 30-cm 
Telescope at the Crow Observatory-Portalegre, and analyse all the photometric and 
radial velocity data published so far. We did not observe significant periodicities in 
the Fourier spectrum of the observed minus calculated (O-C) transit timing and radial 
velocity diagrams (the highest peak having false alarm probabilities equal to 56 per 
cent and 31 per cent respectively) or long term trends. Combining all the available in- 
formation, we conclude that the radial velocity and transit timing techniques exclude 
at 99 per cent confidence limit any perturber more massive than M > 100 M oa rth with 
periods up to ten times the period of the inner planet. We also investigate the possible 
presence of an exomoon on this system and determined that considering the scatter of 
the O-C transit timing residuals a coplanar exomoon would likely produce detectable 
transits. This hypothesis is however apparently ruled out by observations conducted 
by other researchers. In case the orbit of the moon is not coplanar the accuracy of our 
transit timing and transit duration measurements prevents any significant statement. 
Interestingly, on the basis of our reanalysis of SOPHIE data we noted that WASP-3 
passed from a less active (logR hk = —4.95) to a more active (logR hk = —4.8) state 
during the 3 yr monitoring period spanned by the observations. Despite no clear spot 
crossing has been reported for this system, this analysis claims for a more intensive 
monitoring of the activity level of this star in order to understand its impact on pho- 
tometric and radial velocity measurements. 

Key words: techniques: photometric, radial velocities - planets and satellites: indi- 
vidual: WASP-3b - stars: activity 



1 INTRODUCTION 

The field of exoplanets is blessed in these years by an impres- 
sive flow of new exciting discoveries. As the sample of known 
exoplanetary systems increases several interesting character- 
istics become evident posing new challanging issues for the- 
ories of planet formation and evolution. Given their short 
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periods and their large masses the so called Hot-Jupiters 
were the first class of exoplanets being discovered around 
solar type stars (Mayor & Queloz 1995). Those among them 
later found to transit in front of the disk of the star (Char- 
bonneau et al. 2000) also gained a special importance given 
that they allow us to acquire physical informations like the 
radius and the density of the planet which would remain oth- 
erwise inaccessible. At the time of writing this paper, there 
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were 187 known and confirmed transiting planetsQ and 87 of 
them have a period smaller than 10 days and a mass larger 
than 0.7 M jup . 

Several hypothesis were endeavored to explain the ori- 
gin of these objects involving scenarios where these giant 
planets, while originally forming in remote regions of the 
planetary system, were then moved to their actual posi- 
tion by means of different possible mechanisms which can 
be essentially grouped in three broad classes: (i) planet- 
protoplanetary disk interactions leading to inward migra- 
tion of the giant planet (Goldreich & Tremaine 1980; Nelson 
et al. 2000); (ii) planet-planet scattering in mult i- planetary 
systems (Weidenschilling & Marzari 1996; Chatterjee et al. 
2008; Juric & Tremaine 2008); (iii) Kozai-induced migration 
in inclined planetary or binary stellar systems (Kozai 1962; 
Wu & Murray 2003; Fabrycky & Tremaine 2007). While 
the first class of mechanisms would produce in principle 
a smooth migration leading to circularized orbits preserv- 
ing the original alignment between the spin axis of the star 
and the orbital angular momentum axis of the planet, the 
other two processes may result in final eccentric orbits and 
largely misalinged spin-orbit angles. As it was evidenced 
by exploiting the Rossiter-McLaughlin effect (Rossiter 1924; 
McLaughlin 1924, hereafter RM effect), most transiting exo- 
planets have spin-orbit angles perfectly consistent with zero, 
but some of them present surprisingly large misaligned an- 
gles (Triaud et al. 2010). These differences highlight the fact 
that probably all these processes are playing a role in shap- 
ing the structure of these systems (Nagasawa et al. 2008) . 

To further clarify the relative importance of these differ- 
ent theoretical scenarios and understand under which situa- 
tions one mechanism may prevail over the others, some addi- 
tional and important related questions need to be carefully 
examinated. One of them concerns the need to understand if 
these objects are actually isolated or if other planets or even 
stellar companions are gravitationally bounded with the sys- 
tem. Despite the importance of this topic in the framework 
of our understanding of Hot- Jupiter planets, our knowledge 
is still far from being complete. 

In this paper, while attempting to shed new light on 
this problem, we considered the case of the transiting Hot- 
Jupiter WASP-3b, collecting all the photometric and radial 
velocity data acquired so far as well as presenting our new 
photometric measurements. We used this database to in- 
vestigate the presence of an additional companion in this 
system. 

WASP-3b is a Hot- Jupiter planet with a mass of (2.00 ± 
0.09) Mj U p revolving around a main sequence star of spec- 
tral type F7-8V with a period of ~ 1.8 days. Its discov- 
ery was announced in 2008 by the WASP Consortium (Pol- 
lacco et al. 2008) as a result of a photometric campaign 
conducted with the robotically controlled WASP-North Ob- 
servatory located in La Palma and subsequent radial ve- 
locity follow-up obtained with the SOPHIE spectrograph 
at the Observatory de Haute-Provence. The first photom- 
etry of WASP-3b was presented in the discovery paper of 
Pollacco et al. (2008) which used SuperWASP-N together 
with IAC80cm and Keele 80cm telescopes data to refine the 
properties of the transiting object. Two additional transits 
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Table 1. Informations on our observing runs. 



Date 


Epoch 


Texp 


Airmass range 


N. images 


15/05/2009 


196 


90 


2.016-1.002 


178 


13/04/2011 


574 


90 


1.873-1.004 


104 


26/04/2011 


581 


90 


2.055-1.005 


123 


02/06/2011 


601 


150 


1.532-1.024 


96 


20/07/2011 


627 


150 


1.029-2.019 


98 


13/08/2011 


640 


150 


1.032-1.635 


77 


26/08/2011 


647 


150 


1.002-1.644 


77 


08/09/2011 


654 


150 


1.002-1.333 


61 



of WASP-3b were observed by Gibson et al. (2008) with the 
RISE instrument mounted on the fully robotic 2m Liver- 
pool Telescope. Tripathi et al. (2010) observed six transits 
of WASP-3b at the 1.2m FLOW telescope and at the 2.2m 
University of Hawaii Telescope. Joining their results with 
those of Pollacco et al. (2008) and Gibson et al. (2008) they 
concluded that a linear fit to the observed ephemerides was 
not satisfactory, and that either the errors were underesti- 
mated or there was a genuine period variation. Maciejewski 
et al. (2010) presented six new transits gathered at two lm 
class Telescopes (Jena and Rozhen) pointing out that a pe- 
riodic signal was present in the observed minus calculated 
(hereafter O-C) transit timing diagram, and that an outer 
perturbing planet in the system could have best explained 
the observations. Later on Christiansen et al. (2011) dis- 
cussed eight new transits of WASP-3b observed during the 
NASA EPOXI Mission of Opportunity. Despite the high pre- 
cision of their transit timing measurements, these data have 
never been used so far to analyze transit timing variations 
of WASP-3b. Recently Littlefield (2011) reported five ad- 
ditional transit measurements of WASP-3b which were ob- 
served with the 11-inch Schmidt-Cassegrain telescope at Jor- 
dan Hall on the University of Notre Dame Campus. Despite 
the larger uncertainties with respect to previous studies the 
analysis of Littlefield apparently provided an initial mod- 
est support to the hypothesis of Maciejewski et al. (2010). 
Very recently Sada et al. (2012) obtained three additional 
lightcurves of WASP-3b observing with the KPNO visitor 
center 0.5m Telescope and one with the 2.1m KPNO Tele- 
scope. 

Here we present a study of eight new homogeneously 
observed transits of WASP-3b. This paper is structured as 
follows: in Sect. [3] we present our observations of WASP-3b; 
in Sect. |3] we derive the stellar paramters of the host star; 
in Sect.|3j we describe the reduction process; In Sect. [5] we 
describe our analysis of the photometric data; in Sect. [5] we 
present the radial velocity data. In Sect. [7J we describe our 
analysis of the radial velocity data. In Sect. [8] we discuss the 
O-C trasit timing diagram while in Sect. [5] the O-C radial 
velocity diagram. In Sect llOl we discuss our results. Finally 
in Sect. 1 1 1 1 we summarize our results and conclude. 



2 OBSERVATIONS 

The data described here were acquired at the Crow 
Observatory-Portalegre in Portugal. Eight different transits 
of WASP-3b were observed as documented in Table [TJ The 
telescope is a 30 cm aperture Meade LX200 F10, reduced at 
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Table 2. Technical specifications of the acquisition camera. 



Camera Sbig ST8XME specifications 



CCD 


Kodak KAF-1603ME + TI TC-237 


Pixel Array 


1530 x 1020 pix 


CCD Size 


13.8 x 9.2 mm 


Total Pixels 


1.6 million 


Pixel Size 


9/im x 9/mi 


Full Well Capacity 


~100,000 ~e 


Dark Current 


1 - e/pix/sec at °C. 


Readout Specifications 




Shutter 


Elect ro mechanic al 


Exposure 


0.12 to 3600 sec 


resolution 


10 msec 


A/D Converter 


16 bits 


A/D Gain 


2.17 "e/ADU 


Read Noise 


15 ~e RMS 


Binning Modes 


lxl, 2x2, 3x3 


Full Frame Download 


~4 sec 



F5.56 (1668 mm) focal lenght yielding a total field of view of 
~28'xl9. The images were acquired with a Sbig ST8XME 
camera which technical characteristics are reported in Ta- 
ble [2] The pixel scale is 1.1 /pix. The exposure time was 
fixed either to 90 sec or to 150 sec, the overhead was 4 sec 
due to the full frame download time. A total number equal 
to 814 images were acquired and analyzed. All of them were 
in the /-band filter. 



3 DATA REDUCTION 

Bias subtraction and flat fielding were performed with our 
own software in a standard manner. We construct a master 
dark image to identify defective pixels in the image and ap- 
plied a bad pixel correction algorithm which interpolated the 
values of the bad pixels with those of the surrounding pix- 
els. Then we used daophot (Stetson 1987) to derive initial 
aperture photometry and calculate the point spread func- 
tion (PSF) of our images. ALLSTAR was used to refine mag- 
nitude estimates and centroid positions. We then selected 
our best seeing image as astrometric reference frame. Coor- 
dinate transformations among all the other frames and the 
reference were calculated using daomatch and daomaster. 
We took the first ten best seeing images to construct a mas- 
ter high S/N reference frame with MONTAGE2, and a master 
list of objects. After that centroid positions and magnitudes 
were further refined using ALLFRAME (Stetson 1994). Finally 
we rederive aperture photometry for each source after sub- 
tracting the PSF of all the other objects in our images. After 
some experiments we decided to set the aperture radius for 
each frame to 2.3 times the value of the FWHM of the cor- 
responding PSF. 



two distinct normalizations before and after the meridian 
crossing in order to match the photometric zero points. 



3.2 Time stamps 

We report all the mid-exposure times of our measurements 
to the Barycentric Julian Date (BJD) reference frame and 
barycentric dynamical time standard (TDB) using the on- 
line converter provided by Jason Eastman (Eastman, 
Siverd & Gaudi 2010). 



4 STELLAR PARAMETERS 

We used a combined spectrum of WASP-3 to derive spectro- 
scopic stellar atmospheric parameters, including its effective 
temperature and metallicity. The spectrum used is a stacked 
of 8 individual spectra obtained between July and August 
2007 (Pollacco et al. 2008). The spectra were downloaded 
from the OHP-SOPHIE archive. All spectra were obtained 
in the HE mode (R~40 000) placing the fiber B on the sky. 
We used the spectrum in fiber B to subtract any contamina- 
tion light in fiber A (pointing to WASP-3), after correcting 
for the relative efficiency of the two fibers. The final spec- 
trum has a S/N of the order of 100 in the 6500 A region. 

We used the metholology and line-list described in San- 
tos et al. (2004). In brief, after the measurement of the line 
equivalent widths (EWs), the parameters are obtained mak- 
ing use of a line-list of 22 Fel and 9 Fell lines and forcing 
both excitation and ionization equilibrium. We refer to San- 
tos et al. for details. The analysis was done in LTR using a 
grid of Kurucz (1993) model atmospheres and a recent ver- 
sion of the radiative transfer code MOOG Sneden (1973). 
The EWs were derived manually using the IRAF splot task. 

The final obtained stellar parameters are as 
follows: T off =6448±123K, logg=4.49±0.08 dex, 

£i=2.01±0.40kms _1 , and [Fe/H]=-0.02±0.08 dex. As a 
double check, we also independently derived the effective 
temperature of the star using the line-ratio procedure 
described in Sousa et al. (2010). This procedure uses a 
different line-list, and the EWs are measured automatically 
using the ARES (Sousa et al. 2007) code. The effective 
temperature derived using this method is 6432±94 K, in 
perfect agreement with the value mentioned above. These 
values are in agreement with the ones presented in the 
planet announcement paper (Pollacco et al. 2008), who 
derived a temperature of 6400±100K, a surface gravity of 
4.25±0.05 dex, and a metallicity of 0.00±0.20 dex. 



3.1 Corrected lightcurves 

We then constructed the flux ratios between our target 
source (WASP-3b) and several other surrounding compari- 
son stars. We used the first twenty brightest stars in our field 
of view, and calculated a robust weighted average of their 
fluxes after removing any linear differential extinction trend 
as measured in the out-of-transit segments of the lightcurves. 
Since the telescope has a german mount once it crosses the 
meridian it flips around the field of view of 180 degrees. 
Once this event happend we found it was necessary to apply 



5 TRANSIT ANALYSIS 

We modeled the observed transits considering the analytical 
formula of Mandel & Agol (2002). We adopted in particular 
the following parametrization for the planet distance to the 
stellar center normalized to the stellar radius (z): 



2 http: / /astroutils. astronomy.ohio-state.edu/time/utc2bjd. html 
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where G is the gravitational constant, P is the orbital period 
of the planet, p* is the mean stellar density, To is the time of 
transit minimum, Td is the total transit duration (from the 
first to the fourth contact) and r is the ratio of the planetary 
radius to the stellar radius. We fit each lightcurve with the 
Levemberg-Marquardt algorithm (Press 1992). We made use 
of the partial derivatives of the flux loss calculated by Pal 
(2008) as a function of the radius ratio r and the normalized 
distance z. 

The flux F of the star at each given instant of time t 
during the transit (corresponding to the normalized distance 
z) was assumed to be 



F{z(t)) = F MA (z(t))+Ax Air(t)+B, 



(2) 



where Fma is the flux loss predicted by the Mandel & Agol 
(2002) formula, Air(t) is the airmass at the instant t, and A 
and B are two parameters to account for a residual photo- 
metric trend with airmass and a constant zero point offset. 
We therefore assumed six free parameters: the time of tran- 
sit minimum (To), the airmass coefficient (A), the constant 
zero point B, the planet to star radius ratio (r), the transit 
duration (Td) and the mean stellar density (p*). 

We assumed a quadratic limb darkening law which coef- 
ficients were fixed interpolating the tables of Claret & Bloe- 
men (2011) in correpondence to the spectroscopic parame- 
ters of the star. This procedure yielded the following coef- 
ficients in the /-band: gi = 0.2150 for the linear term and 
gi = 0.3034 for the quadratic term. The orbital period of the 
planet was fixed as well at the value of P — 1.846834 days 
(Pollacco et al. 2008). For each iteration the Levemberg- 
Marquardt algorithm calculated the reduced Xred of the fit 
defined as: 



Xred 



\ 



N - 



(Oi - Fi] 



N - N fr 



free 



(3) 



where Oi is the observed flux corresponding to the i-th 
measurement, Fi is the model calculated flux as described 
above, N is the total number of measurements, and Nf ree 
the number of free parameters. The Levemberg-Marquardt 
algorithm found the best solution by means of Xred mini- 
mization. This solution is however only a formal solution, 
the best parameters and their uncertainties were then found 
using a Markov Chain algorithm as described below. 



5.1 Uncertainties of the observations 

To each measurement in our datasets we associated an un- 
certainty accordingly to the photon noise and the read out 
noise as determined by daophot. These errors were then 
added in quadrature to the scatter of the residual fluxes 
of the comparison stars around our derived mean averaged 
values (see Sect. 13. ip . and then rescaled in such a way that 
the best models we fit to the data produced a x 2 ~ 1- As 
pointed out by Pont et al. (2006) the presence of corre- 
lated noise in the data strongly limits the precision of the 



observations. The uncertainties calculated by daophot al- 
ready accounted for some obvious noise correlations. This 
is evident in Fig. [TJ where we plot the uncertainty of the 
measurements as a function of airmass and seeing. This 
ensured that the transits are fitted giving more weight to 
those measurements acquired under the best observing con- 
ditions. Nonetheless, it may well be that the noise in our 
data is correlated also with some other non trivial variables 
that our reduction did not take into account. In order to 
verify this hypothesis we created some mock lightcurves of 
our model-subtracted lightcurves assuming that each simu- 
lated point was distributed normally around zero but with 
a time-dependent dispersion equal to the uncertainty of the 
correspondent real data. We then compared the RMS of the 
real and simulated lightcurves averaged over timescales com- 
prised in between 10 min to 30 min. The average ratio of 
the dispersions of the real to the simulated data (a r /a a ) 
was always smaller than one with the exception of epoch 
574 (cr r /(T s — 1.09), and epoch 627 (a r /a s — 1.04) observa- 
tions. For these two nights we expanded our uncertainties 
by these factors, whereas for the remaining nights we didn't 
apply any other correction. 



5.2 Markov Chain Monte Carlo analysis 

A commonly used approach to derive parameter uncertain- 
ties in exoplanetary literature (e. g. Gazak, J. et al. 2012) 
is based on Markov Chain Monte Carlo analysis. We im- 
plemented our own version of the Markov Chain algorithm 
along the following lines. For each transit lightcurve we cre- 
ated five chains of 10 steps. Each chain is started from a 
point 5-<t away (in one randomly selected free parameter) 
from the best-fitting solution obtained by the Levenberg- 
Marquardt algorithm, where the <j values of the parameters 
considered are those that are obtained by the same transit 
fitting algorithm as described in Sect. [5] The Xoid °f the fit 
of this initial solution is recorded and compared with the 
Xnow obtained in the following step. The following step is 
obtained jumping from the initial position to another one 
in the multidimensional parameter space randomly select- 
ing one of the free parameters and changing its value by an 
arbitrary amount which is dependent on a jump constant 
and the uncertainty a of the parameter itself. Steps are ac- 
cepted or rejected accordingly to the Metropolis-Hastings 
criterium. If Xnow is lower than Xoid the step is executed, 
otherwise the execution probability is P = e _Ax ^ 2 where 
Ax 2 = Xnow — Xoid- I n this latter situation a random num- 
ber between and 1 is drawn from a uniform probability 
distribution. If this number is lower than P then the step 
is executed, otherwise the step is rejected and the previous 
step is repeated instead in the chain. In any case the value 
of the x 2 °f the last step is recorded and compared with 
the one of the following step up to the end of the chain. 
We adjusted the jump constants (one for each parameter) 
in such a way that the step acceptance rate for all the pa- 
rameters was around 25 per cent. The convergence among 
the five separate chains was checked comparing the variances 
within and between the different chains by means of the Gel- 
man & Rubin (1992) statistic. In all cases the values of the 
Gelman-Rubin statistic was within a few percent from unity 
indicating that the chains were converged and well mixed. 
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We then excluded the first 20 per cent steps of each chain 
to avoid the initial burn-in phase, and for each parameter 
we merged the remaining part of the chains together. Then 
we derived the mode of the resulting distributions, and the 
68.3 per cent confidence limits defined by the 15.85th and 
the 84.15th percentiles in the cumulative distributions. 

We run two separate groups of chains first considering as 
free parameters To, A and B while retaining the others fixed 
at the best values obtained by the Levemberg-Marquardt al- 
gorithm. In a second step we instead perturbed p*, Td and r 
while retaining the remaining parameters fixed at the values 
obtained in the first step. Perturbing all the parameters to- 
gether lead in general to unstable convergence in particular 
for the airmass and zero point coefficients so we decided to 
split the procedure in two steps. 



5.3 Mean stellar density 

The mean stellar density is one of the most important 
parameters that can be extracted from transiting planet 
lightcurves (Sozzetti et al. 2007). It is interesting to com- 
pare the stellar density derived from the analysis of our 
lightcurves to the value obtained from the analysis of other 
datasets. 

We then considered the precise transit lightcurve of 
WASP-3b obtained by Tripathi et al. (2010) in the Sloan 
z'-band filter with the University of Hawaii 2.2m Telescope. 
We chose this transit because all the other transits of WASP- 
3b published so far (both by Tripathi and by other authors) 
have been observed with smaller telescopes, and because be- 
ing the observations carried out in a near infrared filter the 
impact of limb-darkening on the transit shape should be 
lower than at shorter wavelenghts. We performed on this 
lightcurve the fit and the statistical analysis previously de- 
scribed obtaining final values for the parameters consistent 
with those of Tripathi et al. (2010). For the mean stellar den- 
sity we obtained p* = 0.50lo'Jg gem -3 . On the contrary the 
analysis of our lightcurves favours a larger density equal to 
= (0.80 ± 0.07) gem -3 by taking the mean average of 
the results reported in Table [3] We note that while these 
estimates are consistent within 2a the value obtained from 
our lightcurves is more similar to the one reported by Pol- 

07 S cm " 



lacco et al. (2008), p* = (0.55jo 5 15 ) P© = 0.77to.o 1 ~~~~ 3 



and Miller et al. (2010) p* = (0.67^ 6 05 ) p = 0.94lg 



+0.07 



gem 



5.4 Results 



In Table [3] we reported for each transit our measured transit 
durations (Td) and planet to stellar radius ratios (r). In 
particular we observed a weighted average transit duration 
Td equal to Td = (158 ± 1) min, which is closer to the value 
reported by Pollacco et al. (2008) T d = (159.8^2 j) min and 
Maciejewski et al. (2010) T d = (161.2 ± 2.3) min, than those 
reported by Tripathi et al. (2010), T d = (168.8 ± 0.7) min 
and Gibson et al. (2008) T d = (165.2±i;jj) min. 

Our weighted average planet to stellar radius ratio 
is equal to r = (0.1061 ± 0.0007) and it is consistent 
with the value reported by Maciejewski et al. (2010) r = 
(0.108 ± 0.003) and Tripathi et al. (2010) for the z filter 
r = 0.1099^0 ooio > but it is larger than the values given 
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Figure 1. Noise correlation with atmospheric indicators relative 
to April 13, 2011. 



by Pollacco et al. (2008) r = 0.1030lHoi° and Gibson et 
al. (2008)r = 0.1014±™. 

Table [4] lists the collection of transit timings of WASP- 
3b presented in published papers, along with our new mea- 
surements. Our timing errors are comprised between 80 sec 
and 233 sec. Note that Maciejewski et al. (2010) transit tim- 
ings are expressed in BJD based on TT (Terrestrial Time). 
The difference with respect to BJD based on TDB is how- 
ever negligible for our purposes (Eastman, Siverd & Gaudi 
2010). Transit timings of Pollacco et al. (2008), Tripathi et 
al. (2010) and Gibson et al. (2008) have been corrected to 
account for the conversion between HJD and BJDtd^E 
We considered all the transit ephemerides presented in 
Table [4] and recalculated the transit period by fitting a 
weighted linear least square model to all the data obtain- 
ing: 



T C (E) = (2454605.5601 ± 0.0002) + 
+ E x (1.846834 ± 0.000001) 



(4) 



where E is the transit epoch. We then subtracted the model 
from the observed ephemerides which gave the (O-C) resid- 
uals presented in Table [4] and in Fig. [3] Our measurements 
are consistent with the calculated ephemerides with the ex- 
ception of those relative to epochs 196 and 574. The reduced 



chi-squared value of the fit (yxP) is equal to 2.30, obtained 
from all the 40 measurements reported in Table U and con- 
sidering 2 degrees of freedom. 

In order to further check our ephemerides measurements we 
applied the barycentric method (Szabo et al. 2006, Oshagh 
et al. 2012). This technique calculates the transit center 
(Tc) as the flux weighted average epoch across the tran- 
sit: 



3 Jason Eastman's Barycentric Julian Date Converter, 
http: / /astroutils. astronomy.ohio-state.edu/time/ 
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Figure 2. Observed transits of WASP-3b along with our best-fitting models and residuals. 
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where /; is the normalized flux, tj the time and N the total 
number of measurements within one transit. Transit fitting 
and barycentric method results agree within 1.9-ct (where 
the a is the average of the uncertainties reported in Ta- 
ble [4} for all the transits with the exception of the transit 
occurring at epoch 640 (3.7-er). In any case, we notice that 
this transit is almost a partial transit which is at the limit 
of applicability of the barycentric method. 



6 RADIAL VELOCITIES 

Radial velocities can be used together with transit timing 
variations to place more stringent constrains on the pres- 
ence of a perturbing object in the WASP-3b system. We 
therefore gathered all the radial velocity measurements of 
WASP-3b publically available from the Exoplanet Orbit 
databas^B These measurements were presented in Pollacco 
et al. (2008), Simpson et al. (2010) and Tripathi et al. (2010) 
and in the following we first introduced in more detail these 
datasets. 



400 



200 - 



o 



-200 



400 



4500 5000 5500 

(BJD TDB -2450000) 

Figure 3. Observed minus calculated times of transit minimum. 
Filled circles denote previous literature results, open circles the 
new measurements presented in this work. 



6.1 Available data sets 

Pollacco et al (2008) obtained seven radial velocity mea- 
surements of WASP-3 using the SOPHIE spectrograph at 
the 1.93-m telescope at Haute-Provence Observatory. The 
observations were performed between 2007 July 2-5 and Au- 
gust 27-30. All the measurements were acquired outside the 
transit. Simpson et al. (2010) acquired 26 spectra of WASP-3 
during the transit occurring on the night of 2008 September 
30. The observations were also obtained with the SOPHIE 



4 |http://e xoplanet s.org/| 
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spectrograph at the 1.93-m telescope at Haute-Provence Ob- 
servatory. These authors also reanalysed the seven measure- 
ments presented in Pollacco et al. (2008) based on an up- 
dated version of the SOPHIE pipeline. We therefore decided 
to use these data in our study and not the original data pre- 
sented in Pollacco et al. (2008). 

Tripathi et al. (2010) obtained 33 radial velocity mea- 
surements of WASP-3b with the High Resolution Spectrom- 
eter (HIRES; Vogt et al. 1994) on the Keck I 10 m telescope 
at the W. M. Keck Observatory on Mauna Kea. The ob- 
servations were acquired both during the transit (on 2008 
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Table 3. Best-fitting transit parameters. 



Date 


Epoch (E) 


Slope (s) 


Constant (c) 




few 


Duration (Tj) 
(min) 


Mean stellar density (/?*) 
(gem -3 ) 
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0.00060i°;°W60 


„ _ _„ „_l_n 00077 
f) QQQl 5T u - UUUf ' 
u - yyyio -0. 00094 
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1.00068l8;22?M 
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Table 4. Collection of transit timing measurements of WASP-3b in chronological order. 



Epoch 


Time of transit minimum 


A-(BJD TDB ) 


A+(BJD TDB ) 


(0 - C) 


A~(0 - C) 


A+(0 - C) 


Reference 




(BJDtdb-2450000) 


(days) 


(days) 


(sec) 


(sec) 


(sec) 




-250. 


4143.85104 


0.00040 


0.00040 


-46. 


35. 


35. 


Pollacco et al. (2008) 


-2. 


4601.86588 


0.00027 


0.00027 


-44. 


23. 


23. 


Tripathi et al. (2010) 


0. 


4605.56030 


0.00035 


0.00035 


21. 


30. 


30. 


Gibson et al. (2008) 


12. 


4627.72172 


0.00031 


0.00031 


-30. 


27. 


27. 


Tripathi et al. (2010) 


18. 


4638.80403 


0.00031 


0.00031 


83. 


27. 


27. 


Tripathi et al. (2010) 


30. 


4660.96509 


0.00021 


0.00021 


1. 


18. 


IS. 


Tripathi et al. (2010) 


40. 


4679.43269 


0.00050 


0.00050 


-63. 


43. 


43. 


Christiansen et. al. (2011) 


41. 


4681.27911 


0.00040 


0.00040 


-99. 


35. 


35. 


Christiansen et. al. (2011) 


42. 


4683.12740 


0.00035 


0.00035 


27. 


30. 


30. 


Christiansen et. al. (2011) 


43. 


4684.97486 


0.00027 


0.00027 


81. 


23. 


23. 


Christiansen et. al. (2011) 


44. 


4686.82053 


0.00059 


0.00059 


-19. 


51. 


51. 


Christiansen et. al. (2011) 


46. 


4690.51381 


0.00055 


0.00055 


-53. 


48. 


48. 


Christiansen et. al. (2011) 


47. 


4692.36117 


0.00043 


0.00043 


-7. 


37. 


37. 


Christiansen et. al. (2011) 


48. 


4694.20711 


0.00042 


0.00042 


-84. 


36. 


36. 


Christiansen et. al. (2011) 


59. 


4714.52284 


0.00036 


0.00036 


-36. 


31. 


31. 


Gibson et al. (2008) 


194. 


4963.84436 


0.00072 


0.00072 


-128. 


62. 


62. 


Tripathi et al. (2010) 


194. 


4963.84563 


0.00055 


0.00055 


-18. 


48. 


48. 


Sada et al. (2012) 


196. 


4967.53651 


0.00057 


0.00085 


-259. 


49. 


73. 


This work 


201. 


4976.77365 


0.00051 


0.00051 


-3. 


44. 


44. 


Tripathi et al. (2010) 


236. 


5041.41271 


0.00049 


0.00049 


-14. 


42. 


42. 


Maciejewski et al. (2010) 


249. 


5065.41995 


0.00059 


0.00059 


-152. 


51. 


51. 


Maciejewski et al. (2010) 


256. 


5078.34873 


0.00058 


0.00058 


-71. 


50. 


50. 


Maciejewski et al. (2010) 


269. 


5102.35933 


0.00056 


0.00056 


81. 


48. 


48. 


Maciejewski et al. (2010) 


289. 


5139.29713 


0.00049 


0.00049 


178. 


42. 


42. 


Maciejewski et al. (2010) 


379. 


5305.51082 


0.00039 


0.00039 


60. 


34. 


34. 


Maciejewski et al. (2010) 


403. 


5349.83457 


0.00039 


0.00039 


37. 


34. 


34. 


Sada et al. (2012) 


403. 


5349.83182 


0.00039 


0.00039 


-200. 


34. 


34. 


Sada et al. (2012) 


404. 


5351.68320 


0.00110 


0.00110 


192. 


95. 


95. 


Littlefleld (2011) 


430. 


5399.69990 


0.00150 


0.00150 


108. 


130. 


130. 


Littlefleld (2011) 


436. 


5410.78020 


0.00130 


0.00130 


47. 


112. 


112. 


Littlefleld (2011) 


450. 


5436.63590 


0.00080 


0.00080 


49. 


69. 


69. 


Littlefleld (2011) 


456. 


5447.71550 


0.00080 


0.00080 


-72. 


69. 


69. 


Littlefleld (2011) 


574. 


5665.64627 


0.00069 


0.00056 


305. 


60. 


48. 


This work 


581. 


5678.57065 


0.00106 


0.00087 


6. 


92. 


75. 


This work 


592. 


5698.88358 


0.00060 


0.00060 


-188. 


52. 


52. 


Sada et al. (2012) 


601. 


5715.50608 


0.00074 


0.00060 


-102. 


64. 


52. 


This work 


627. 


5763.52552 


0.00070 


0.00047 


50. 


60. 


41. 


This work 


640. 


5787.53379 


0.00080 


0.00080 


1. 


69. 


69. 


This work 


647. 


5800.46112 


0.00113 


0.00170 


-43. 


98. 


147. 


This work 


654. 


5813.38792 


0.00098 


0.00080 


-133. 


85. 


69. 


This work 
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Table 5. Adopted values for the p and q parameters and for the 
limb darkening coefficients entering in the Rossiter-McLaughlin 
model. 



p 


q 


91 


92 


Reference 


1.51 


0.44 


0.596 


0.215 


Tripathi et al. (2010) 


1.72 


0.00546 


0.69 


0. 


Simpson et al. (2010) 



June 19, 21 and on 2009 June 3) and outside the transit on 
several other nights in 2008 and 2009. 



7 ANALYSIS OF THE RADIAL VELOCITY 
DATA 

We analyzed simultaneously all the radial velocity measure- 
ments presented both by Tripathi et al. (2010) and Simp- 
son et al. (2010). Since many measurements were acquired 
during the transit, we fit the data with a model describing 
both the Keplerian motion of the host star and the Rossiter- 
McLaughlin anomaly. On one hand we modeled the Keple- 
rian motion as: 



RV = K 



cos u + k 
VI - h 2 - fc f 



+ 7 



(6) 



where K is the radial velocity semi-amplitude without the 
contribution of the eccentricity e, k = e cosoj, h = e sinw, 
7 is the barycentric radial velocity and u = v + ui is the 
true argument of latitude with v the true anomaly and uj 
the argument of the pericenter. 

On the other hand we accounted for the Rossiter- 
McLaughlin (RM) anomaly following Hirano et al. (2010) 
but including an improved treatment of the RM effect dur- 
ing the partial phases of the transit (as detailed in Appendix 
B): 



RV RM = -df x v p 



V - 



(-) 



where df is the flux loss due to the transit of the planet 
in front of the disk of the star, which we modeled as in 
Mandel & Agol (2002), p and q are two parameters related 
to modellization of the Rossiter-McLaughlin effect as pro- 
posed by Hirano et al. (2010) and were fixed to the values 
adopted by Tripathi et al. (2010) and Simpson et al. (2010) 
as reported in Tabled v p is the average velocity of the star 
below the area occulted by the planet (Hirano et al. 2010 
and Appendix B), and vsini is the rotation velocity of the 
star. 

We considered as free parameters: K, h, k, X (the spin- 
orbit angle), usini, and 7. Both Tripathi et al. (2010) and 
Simpson et al. (2010) distinguished two different groups of 
data in their own dataset to account for possible systematic 
radial velocity variations during their observing runs. We 
decided to perform the fit twice, first following the analy- 
sis of those authors and therefore allowing for a total num- 
ber of four different barycentric radial velocities. Then, we 
redid the fit considering only two different barycentric ra- 
dial velocities for the Tripathi et al. (2010) and Simpson et 
al. (2010) datasets. This second approach was intended to 



check for possible long-term variations among the RV resid- 
uals that could have been canceled out by the adoption of 
a larger number of free parameters. Then, in the end, we 
considered either nine or seven free parameters for the two 
fits respectively. We notice that Tripathi et al. (2010) ex- 
cluded from the fit three data points which were presenting a 
clearly deviant radial velocity with respect to the remaining 
measurements. This radial velocity spike was ultimately at- 
tributed by the authors to residual moonlight unexpectedly 
leaking into the spectrograph and therefore we neglected 
them hereafter. Additionally Tripathi et al. (2010) added in 
quadrature to the uncertainties of their data a value equal 
to 14.8 ms -1 to account for jitter noise. Since, however, it 
is not clear which is the origin of this noise, we didn't apply 
this correction. 

The convergence toward the best-fit solution was ob- 
tained by means of a Levenberg-Marquardt algorithm, and 
the uncertainties and the best fit values of the parameters 
by means of a Markov Chain Monte Carlo analysis as done 
for the photometric data. The radial velocity measurements 
after subtraction of the barycentric velocities, along with 
the best fit model and distinguished in the four (two) dif- 
ferent groups of data to which the different 7 were applied 
are shown in Fig. 2] (upper panels). The result after the 
subtraction of the RM anomaly is shown in the middle pan- 
els of Fig. [4] and the residual velocities after subtracting 
the Keplerian orbit also are shown in the bottom panels. 
Our best-fit parameters are given in Table [6] along with 
the values obtained by other authors. Our best-fit model 
corresponds to a y / xf = 1.5 for the four 7 solution and 
to VaI = 1-6 f° r the two 7 solution. Our results are in 
agreement with the literature values. The barycentric ve- 
locities for the case of the four 7 solution are consistent 
with those derived by Tripathi et al. (2010) and Simpson 
et al. (2010), and in the case of the 27 solution our val- 
ues for each dataset are in between the results reported 
by those authors for their own data. We obtained a value 
of the spin-orbit angle consistent with zero. We also notice 
that the rotation velocity we obtained for the four 7 solution 
(vsin(i) = 13.9lo 5) * s perfectly consistent with the result of 
Miller et al. (2010) implying v sin(i) = 13.9±g For the two 
7 solution we obtained instead a larger value of the rotation 



velocity (vsin(i) 



14.5 



+0.3\ 



8 ANALYSIS OF THE (O-C) TRANSIT 
TIMING DIAGRAM 

In a first step, we analyse whether or not a quadratic de- 
parture from a linear fit is present in the transit timings. 
This could result from the direct interaction with a per- 
turber on an extended orbit (Borkovits et al. 2011), or from 
the light travel timing produced by the motion of the star 
also induced by a hypothetical distant companion (Mon- 
talto 2010). This test can be performed following the ap- 
proach of Pringle (1975) which measures the improvement 
of a fit by a quadratic parabola with respect to a simpler 
one by a straight line. However, the quadratic coefficient 
obtained by least square is already zero within the errorbars 
(—0.9 ±3.5) x 10 -9 days. We thus conclude that there is no 
significant long term quadratic trend in the data. 

Then, we follow the analysis of Maciejewski et al. 
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Figure 4. Radial velocity fit once considering four and two values of 7 as described in the text (left and right panels respectively). 
Top panels: radial velocity measurements along with our best fit models. Barycentric radial velocities are subtracted. Middle panels: all 
radial velocity measurements after subtraction of the barycentric velocities and the Rossiter-McLaughlin anomaly. Lower panels: radial 
velocity measurements after subtraction of the barycentric velocities, the Rossiter-McLaughlin anomaly and the Keplerian orbit. 
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Table 6. Best-fit parameters obtained from our reanalysis of the radial velocity measurements (TW=This work), and from the study 
of Simpson et al. (2010, SI10), Tripathi et al. (2010, TRIO), Miller et al. (2010, MHO) and Pollacco et al. (2008, PO08). 



v sin i 
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71 
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(2010). We compute a Lomb-Scargle periodogram on the 
Transit Timing Variation (TTV) signal in order to detect a 
periodic oscillation that would reflect the perturbation of a 
close-in undetected body in the system. For that purpose, 
we use the generalized version (GLS) of the Lomb-Scargle 
periodogram (Zechmeister & Kiirster 2009). Basically the 
GLS fits a sinusoid to the data for each frequency by using 
the least square method, as the Lomb-Scargle algorithm, 
but in addition to that it allows for the presence of an ad- 
ditional constant term. False-alarm probabilities (FAP) are 
estimated by computing GLS periodograms on a large num- 
ber of sets of artificial observations in which, for each epoch 
where a transit has been observed, we replace the measured 
(O-C) by a random value normally distributed around zero 
with a standard deviation equal to the uncertainty of that 
point. The number of periodograms containing a peak with 
an power above a given threshold out of the total number 
of trials represents our estimation of the FAP for that given 
threshold. 

Fig. [5ja) shows the GLS periodogram obtained when 
considering only the transits used by Maciejewski et al. 
(2010, Fig. 3). We obtain a dominant peak at a frequency 
/ttv = 0.0145 cycle P , which corresponds to Pttv = 
127 days and an power of 0.61, as shown by the arrow, in 
complete agreement with the result of Maciejewski et al. 
(2010). Nevertheless the false-alarm probability associated 
to that power is 27%. There is thus more than one chance 
out of 4 for this peak to be fortuitous. For the sake of com- 
pleteness, the false-alarm probability thresholds of 0.1, 10~ 2 , 
and 10~ 3 are represented by three horizontal lines in the two 
panels of Fig. [S] We did again the same analysis with all the 
data of the table [4] The results are displayed in Fig. [5jb) . 
In that case, the peak with the highest power is now at 
/ttv = 0.0201 cycle P" 1 with a FAP equal to 56%. It thus 
seems that the TTV signal does not contain any significant 
periodic oscillations. 



9 ANALYSIS OF THE (O-C) RADIAL 
VELOCITY DIAGRAM 

We initially checked for the presence of either a linear or 
a quadratic term in the (O-C) radial velocity residuals by 
using the Pringle (1975) test. We considered the 2y solution 
and obtained that in both cases the coefficients are consis- 
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Figure 5. Periodograms of the TTV signal, (a) Considering only 
the data used in Maciejewski et al. 2010. (b) Using all the data 
present in Tab. [4] The horizontal lines (from bottom to top) give 
the FAP thresholds 0.1, 10" 2 , and 10" 3 . 





Figure 6. Periodogram of the radial velocity residuals consider- 
ing the 47 solution (left) and the 27 solution (right). The continu- 
ous line denotes the 1% FAP while the dashed line the time-span 
of the observations. 



tent with zero (0.00001 ± 0.00012 for the quadratic term) 
and (-0.030 ± 0.035 for the linear term). 

A GLS periodogram was then computed. We consid- 
ered initially the case of the residuals obtained fitting the 
47 solution (Sect. [7}. As seen in Fig. [6] no significant peaks 
can be found in the periodogram. The highest peak is at 
0.35 days and has a FAP of 31%. Alternative, considering 
the residuals obtained by the 27 solution, we obtained the 
highest peak at 0.36 days with a FAP=39%. The FAP of 
the peaks are estimated with a bootstraping method in the 
same way as in the previous section. The only difference is 
that the artificial data are made by shuffling (with repeti- 
tion) the residuals instead of drawing random values from a 
normal distribution. 
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10 DISCUSSION 

Since the TTV signal does not present any long term vari- 
ations, nor short period oscillations, one cannot assert that 
the system actually contains an additional planet. And, if 
such a companion does exist, its orbital parameters and its 
mass are poorly constrained due to the lack of expected pat- 
terns in the current available data. Nevertheless, the resid- 
uals of the fit are quite large (^/x^ ~ 2.30 and ^ 1.5 
for transit and RV respectively). As noticed also by other 
investigators in the past these values do not indicate a satis- 
factory fit to the (O-C). The model is thus not complete. 
Several causes can be at the origin of this result among 
which we consider here stellar activity, the presence of an 
additional planet or exomoon and underestimated timing 
uncertainties. 



10.1 Stellar activity 

A possible source of TTVs can be the activity of WASP- 
3. Indeed, the existence of spots on the surface of the star, 
partially covered by the planet during transits, should pro- 
duce fluctuations in the luminosity leading to some errors 
in the determination of the times of transit minimum (e. g. 
Sanchis-Ojeda et al. 2011; Oshagh et al. 2012). Moreover, 
the spots, if they exist, should not be the same between the 
beginning and the end of the observations given the long 
time span that has been covered (~ 4.5 years). This would 
explain why no periodic oscillation is detected. Tripathi et 
al. (2010) reported fractional transit depth variations of the 
order of 7%, even if the same authors were not confident 
whether these variations were genuine or due to systematics 
in their data. In addition, they report a mean logR HK =-4.9 
from their Keck spectra taken in 2008-2009. 

On the other hand, we reanalysed the spectra taken 
with SOPHIE. From the 2007 observations (Pollacco et 
al. 2008) we derived a logR HK value of -4.95, whereas the 
2009-2010 observations (Simpson et al. 2010) provided a 
higher value for the activity index of log R HK =-4.80. There- 
fore it appears that the mean activity level of the star 
changed during these years approaching an active phase in 
2010. Once considering also the upper limit on the stellar 
age and the rotation period reported by Miller et al. (2010, 
age< 2 Gyr and P ro t = 4.3 days) the presence of active 
regions on this star may not appear a rare circumstance. 
Despite no clear evidence of starspots crossing has been re- 
ported yet, this analysis clearly claims for a more intensive 
monitoring of the activity level of WASP-3 in order to un- 
derstand its impact on photometric and radial velocity mea- 
surements. 



10.2 Additional planet 

We now give some constraints on the mass of a hypothetical 
planetary perturber. For that, we use both the dispersion 
of the O-C radial velocity and photometric diagrams. For 
what concerns the radial velocity residuals we adopted here 
the results coming from the 27 solution, since no significant 
difference was found adopting instead the 47 solution. 



10.2.1 Radial velocities 

In the literature, two main approaches are used to find the 
detection limits in radial velocity data. One is based on 
X 2 - and F-tests (e.g. Lagrange et al. 2009, Sozzetti et al. 
2009), another is based on a periodogram analysis (Cum- 
ming, Marcy & Butler 1999, Endl et al. 2001, Cumming 
2004, Narayan, Cumming & Lin 2005). Here, the second 
approach was chosen due to the number of measurements 
which is considered high enough for a reliable periodogram 
analysis. 

For each period, a fake eccentric planetary signal is in- 
serted in the data, while the original data is treated as ran- 
dom noise. On these new RV series, the power (in the pe- 
riodogram) is calculated. The semi-amplitude of the fake 
signal is changed untill the FAP level is reached for all ec- 
centricities e, times of periastron T c and longitudes of pe- 
riastron vj. In this paper, a FAP of 1%, determined with 
1000 shuffled time series, is used. Fake signals are tested for 
periods P between 1 and 20 days. The orbital elements of 
the eccentric signals range, in 10 steps, as follows: ^ e ^ 1, 
^ T c ^ P and ^ vo ^ 2n. The final semi-amplitude can 
be transformed in planetary mass and expresses the lower 
limit for detectable planets at that period with these data. 

/ pj\ /2\ 1/3 

Afpsini = 1.2 ■ HT 3 AVl - e 2 ( (8) 

with the planetary mass in Earth mass, the semi-amplitude 
K in m/s, the period P in days, the stellar mass M* in so- 
lar masses and the gravitational constant G in m 3 kg -1 s -2 . 
Therefore the continuous grey line in Fig. [7]denotes the limit 
in the perturber mass beyond which a signal would have 
been detected in the radial velocity data with a confidence 
limit equal to 99%. The dashed lines shows the 1-a uncer- 
tainty range of the radial velocity detection limit. 

10.2.2 Photometry 

We exclude compact systems which lead to unstable evolu- 
tions. Only circular and coplanar systems have been con- 
sidered since they provide the strongest constraints and be- 
cause the projected spin-orbit angle measure on WASP-3b 
by the Rossiter-McLaughlin effect is compatible with zero 
as demonstrated above, suggesting that if a planetary com- 
panion exists, the system is likely coplanar. 

The radial velocity measurements are used to exclude 
any perturbers that would induce RV signal with an ampli- 
tude larger than the 1% FAP threshold. The same exercise 
has been performed with the O-C of the transit timing mea- 
surements. For that, we simulated a large number of O-C 
on a grid of parameters of the perturber. We considered 400 
periods P pcrt ranging between 1.5 and 10 Ptransit (where 
-Ptransit = P is the period of the transiting planet) , and 100 
masses M por t evenly distributed in logarithm between 0.01 
and 5 Mj. For each period and mass, 60 simulations are per- 
formed with different initial longitudes between and 360°. 
Among the 60 simulations, the one giving the periodogram 
with the lowest maximum amplitude is kept. If this ampli- 
tude is above the 1% FAP threshold determined in section [8] 
then the corresponding perturber should have been detected 
in the periodogram of the O-C (see Fig. [5j, otherwise the 
perturber can exist but it is not detectable. 
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Figure 7. Left: constraints on the maximal mass of a perturber: 99% confidence detection limits from periodograms of transit timing 
(black) and radial velocity data (grey). The dashed lines denotes the uncertainty range of the radial velocity detection region. The 
hatched lines indicates unstable orbits. Right: reduced yjxjj. values resulting from the fit of the observed O-C timing residuals with our 
model. The calculation is extended over the radial velocity undetectability region, but we also show the regions that produce a 3<r in 
black. Notice that these regions overlap with the transit timing detectability regions of the left panel. 



Figure [Jj (left) shows the results. The hatched region, 
which extends up to a period ratio of 1.5, delineates the 
chaotic orbits which are excluded. The boundary of this re- 
gion has been derived from the stability criterion of (Glad- 
man 1993). The gray curve fixes the limit of the perturber's 
mass from the radial velocity measurements. Any perturber 
in the gray region would induce a periodic RV signal with 
a significant amplitude. And finally, the black regions de- 
lineate the perturbers that produce TTVs with significant 
oscillating terms. Combining all the information, it turns out 
that the radial velocity technique excludes any perturber as 
massive as Jupiter up to a period ratio of 10, and the TTV 
measurements provide stronger constraints close to mean 
motion resonances. 

In a second step, we focused on perturbers that do 
not produce any significant periodic RV signals nor periodic 
TTV signals. Such perturbers are located in the white re- 
gion of Fig. [JJ (left). For those perturbers, we check whether 
they can reproduce the O-C transit timing diagram or not. 
In that purpose, we use the same grid of parameters as in 
Fig. [7] (left), and for each initial conditions, we compute the 
expected TTVs and the reduced chi-square with respect to 
the observations. The results are displayed in Fig. [7] (right). 
The hached and the gray regions are the same as in Fig. [7] 
(left). The black regions correspond now to simulated TTVs 
above the 3-cr level. This threshold is obtained from the x 2 - 
distribution with 33 degrees of freedom. It corresponds to 
Xr = X?.min(l + 1-821), or y^f? = 2.95. The colour scale 
represents the \/xf from the lowest values in red up to the 
3-cr threshold in dark violet. The red circle shows the best fit 
to the observation with \JXr — 1-76. However, such a per- 
turber, with a mass of 0.63 Mj should have been detected 
in the radial velocity analysis. The best fit within the unde- 
tectable perturbers is just below the RV detection threshold 
with Afport = 0.41Mj and P pcrt = 5.63P t ransit, but the cor- 
responding reduced chi-square is only yxT = 1-83. The im- 
provement is very weak. Moreover, from Fig. [7J (right), one 
can see that such values of the reduced chi-square are spread 
more or less randomly within the undetectable region. 



As noted by Maciejewski et al. (2010), the presence of 
an outer companion less massive than WASP-3b but still on 
a short period orbit would make the system quite unusual. 
Multiplanetary systems containing at least a Jupiter-mass 
planet are indeed much wider, and the less massive planet 
is usually the closest to the star (e. g. Lissauer et al. 2011). 



10.3 Exomoon 

An exomoon is also supposed to generate a periodic oscil- 
lation in the TTV. However, Maciejewski et al. (2010) have 
already discarded this hypothesis since the transits do not 
show any duration variations shifted in phase by 7r/2 with 
respect to the timing variations. Here, we perform a more 
detailled analysis based on the results of Kipping (2009). 

First of all, we check that an exomoon can have a stable 
orbit. If the moon is less than twice as dense as the planet the 
minimum distance of the moon is set by the Roche limit. Let 
£ be the semi-major axis of a hypothetical satellite divided 
by the Hill Radius, i.e. a s = £Rh where a s is the semi- 
major axis of the moon and Rh the Hill Radius. According 
to Kipping (2009), £ should satisfy the following inequality 

Xmin < i < ~ , (9) 

where £ min = 1/186 x (M 3 /M®)-°' 063 (P/ld)- 2/3 represents 
the Roche limit. In this expression, M s is the mass of the 
satellite and P is the orbital period of the planet. For an exo- 
moon of the mass of the Earth's Moon, we get f m i n = 0.0047, 
and for an exomoon of the mass of the Earth, £ m i n = 0.0036. 
In both cases Xmin is lower than 1/3. For moons which den- 
sity is more than twice that of the planet the Roche limit 
would be inside the planet, therefore the minimum distance 
would correspond to the planetary radius and the above in- 
equality would be automatically satisfied. Therefore an ex- 
omoon can exist on a stable orbit around WASP-3. 

Then, we estimate the maximal RMS amplitude (<5ttv) 
of TTV that an exomoon on a coplanar circular orbit can 
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produce. According to Kipping (2009), this amplitude is 
given by 

^^TMm)"^--^' (10) 

with fi = M s /Mp S and M ps is the sum of the planet mass and 
the satellite mass. Without any constraint on /x, the maxi- 
mum of the product fi(l — /i) 1 ^ 3 is attained for = 3/4, and 
is equal to 3/4 1 / 3 . However, by definition, the Moon should 
have a lower mass than the planet. If the transit lightcurves 
are those of the planet, \i should be lower than 1/2, and 
probably much lower. But let us assume that fj, — 1/2, 
this will provide the upper limit of 5ttv- If the planet 
has an exomoon, the Keplerian orbit derived in the previ- 
ous section is that of the planet-satellite barycenter around 
the star. Thus the fitted mass corresponds to M sp . Using 
M 3p /AL = 1.5 x 10~ 3 , one obtains 

For £ = £max = 1/3, this leads to <5ttv 3.1 min. The 
result is larger than the observed RMS of the O-C which 
is equal to 1.1 min. Thus, a "satellite" as massive as the 
planet (fi = 1 /2) is able to produce significant TTV with an 
amplitude comparable to the observed one. 

We now assume that the observed O-C is only due 
to an hypothetical satellite. Expecting that this satellite 
should have a much smaller mass than the planet, we de- 
rive its mass for £ = £ ma x = 1/3 such that <5ttv = 1.1 min. 
One gets /j = 0.14, which corresponds to a mass ratio of 
M a /Mp = 0.17, or M 3 w 0.35Mj. Thus, the lowest massive 
satellite, on circular orbit, that can account for the observed 
RMS of the O-C is still large. Assuming the same density 
as a giant planet like Jupiter, the radius of this satellite 
should be 0.46-Rj. Such a big satellite, if it existed, would 
produce detectable transits in the lightcurve. We also no- 
ticed that a search for additional transiting objects in the 
NASA EPOXI mission WASP-3 lightcurve resulted in a 
null detection (Ballard et al. 2011). However, we note that 
accordingly to Domingos et al. (2006) the value of the crit- 
ical semi-major axis could approach the value of 0.9309 in 
units of the Hill radius, for retrogade moons. In this case 
we can obtain a more stringent constraint on the maximum 
mass of the moon which could be equal to M 3 /M p = 0.04 
or M 3 « 0.07 Mj. In Fig. [8l we also show our transit du- 
rations against the O-C residuals. In case of an exomoon 
being responsible for the claimed TTVs, we should expect 
the observations to trace an ellipse in this diagram since 
the TTVs and the TDVs (Transit Duration Variations) pro- 
duced by an exomoon are shifted in phase by 7r/2 (Kipping 
2009). For illustration, we overplot the expected signal that 
the above-mentioned prograde satellite (M s — 0.35A/j and 
X = 1/3) should generate. Evidently given the large error- 
bars of our measurements it is not possible to explore this 
possiblity, and additional more accurate measurements are 
required for this analysis. 

10.4 Underestimated uncertainties 

It is difficult to ascertain up to which level different instru- 
ments, observing conditions, reduction and transit fitting 
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Figure 8. Transit durations against O-C transit timing residuals 
for the new transits presented in this work. Overplot a represen- 
tative signal produced by an exomoon having M s = 0.35Mj and 
period x = 1/3. 



procedures may affect the results reported in Table [4] To 
address this point it would be necessary to homogeneously 
reduce and analyze all the data collected so far by all the 
different groups, an approach which is not easy to put into 
practice. In principle all the transits considered here were 
presented in referred journals and this ensures that accurate 
procedures like those ones reported here have been applied 
to estimate transit timing errors. It is our opinion however 
that error underestimation cannot be completely ruled out. 
New observations will be certainly welcome to clarify this 
problem. 



11 CONCLUSIONS 

In this work we provided a throughout analysis on the pres- 
ence of additional bodies in the WASP-3 system. This anal- 
ysis serves to improve our understanding of close-in Jupiters 
and in particular to clarify if these planets are indeed iso- 
lated or not. 

In addition to present eight new transits of WASP-3b 
acquired at the Crow-Observatory-Portalegre in Portugal we 
reanalized all the photometric and radial velcocity measure- 
ments aquired so far for this system. We concluded that 
there is no convincing evidence of additional planetary com- 
panions in this system; both the transit timing and the ra- 
dial velocity residuals do not present significant periodicities 
(FAP=56% and FAP=31% for transit and radial velocity in 
the best case scenario respectively) nor long term trends. 

Combining all transit timing and radial velocity infor- 
mations, we obtained that any perturber more massive than 
M > lOOM® and with period up to ten times the period of 
the inner planets is excluded at 99% confidence limit. 

We also investigated the possible presence of an exo- 
moon on this system and determined that considering the 
scatter of the O-C transit timing residuals a coplanar exo- 
moon would likely produce detectable transits, an hypothe- 
sis that can be ruled out by observations conducted by other 
researchers. In case the orbit of the moon is not complanar 
the current accuracy of transit timing and transit duration 
measurements prevents to make any significant statement. 
For retrogade moons the maximum mass allowed at the crit- 
ical semi-major axis is around 0.1 Mj. 
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Finally on the basis of our reanalysis of SOPHIE data 
we noted that WASP-3 passed from a less active ( log R hk = 
—4.95) to a more active ( log R hk = —4.8) state between 
2007 and 2010. Despite no clear spot crossing has been re- 
ported for this system so far we therefore pointed out the 
need for a more intensive monitoring of the activity level of 
this star in order to understand its impact on photometric 
and radial velocity measurements. 

Our lightcurves are made available through the on-line 
version of this journal. 
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Figure 9. Definition of the ipQ, (f>Q, e and f angles introduced in 
the text. 1 



of the planet. Eliminating in the first equation above the 
impact parameter derived from the second, using the third 
Kepler law and introducing the definition of the mean stellar 
density 



APPENDIX A: THE NORMALIZED PLANET 
DISTANCE 

Here we derive the expression for the normalized planet dis- 
tance z presented in Eq. 1. Assuming a circular orbit and 
a constant projected velocity (v) of the transiting planet on 
the plane of the sky, in any given instant t during the transit 
the normalized distance z can be written as 



we obtain Eq. 1: 



z\t) = 



3P 



2/3 



(f - To) 2 - 



+ 



+ (1 + r) 2 



+ 



t-T 



where b is the impact parameter, R± is the radius of the star 
and To is the time of transit minimum. At the time of the 
first or the fourth contact we have 



(r + 1) 2 



T d 
2R t 



where r is the ratio of the radius of the planet to the radius 
of the star, and Td is the total transit duration (from the 
first to the fourth contact). Assuming the projected velocity 
during the transit to be identical to the orbital velocity we 
can write 



where G is the gravitational constant, M* is the mass of the 
star, a is the semi-major axis and we neglected the mass 



APPENDIX B: RM EFFECT DURING THE 
INGRESS AND THE EGRESS 

In this Appendix we introduce a new analytic representation 
of the Rossiter-McLaughlin effect valid during the ingress 
and the egress of the transit, that is between the first to 
the second contact and between the third to the fourth con- 
tact. During these phases the formula presented by Hirano 
et al. (2010) accounts for the velocity of the star below the 
disk of the planet considering the value of the velocity at 
the center of the disk of the planet, and it is therefore valid 
for the small planets approximation. We instead integrated 
the velocity profile below the disk and calculated the aver- 
age velocity which makes our approach consistent with the 
calculation of Hirano et al. (2010) for the remaining phases 
of the transit. This derivation is based on the method de- 
scribed in (Pal 2012). 

Therefore if we define X and Y as the coordinates of the 
center of the planet at a given instant during the transit, 
accordingly to the choice of parameters we adopted in this 
paper we have 
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v _ (1-esina;) ( 8n 2 G y /3 



3P 



1/3 (t - To) , 



87r2Gy /3 n 2/3^r d x 2 
3P 



During the full transit pahse, Eq. 10 reduces to 

v v 



(13) 



which is equal to Eq. A8 of Hirano et al. (2010). 



Then if A is the spin-orbit angle projected on the plane of 
the sky, in the rotated coordinate system which vertical axis 
is aligned with the projected spin axis of the star, the coor- 
dinates x and y are given by 

x — X cos A — Y sin A , 
y — X sin A + Y cos A . 



Let ipo, 4"i, e and £ be defined as in Fig. ?? with respect to 
the xy rotated coordinate system, then 

■00 = atan 2 (y,x) , 

Assuming that the radius of the star is normalized to unity, 
one gets 

1 -r 2 +z 2 



2z 



£ = acos 



2rz 



where 

2 = \J x 2 + y 2 . 

Denning now the angles ip a , ipb, 4>a an d 4>b as 

ipa = ipo - e , 
i>b = i>o + e , 

<t>a = <f>0 - £ , 

4>b = 00 + £ , 

and defining the following functions 

1 2 1 2 ( 1 

A(a, /3, 7) = -q /3sin7H — a/3 I 7 -I — sin 27 



+ \P 3 (7 +i sin 3 2 7 



S(a,/3,7) = -a^sin7+-/3 I 7+ -sin 27) , 



the subplanet velocity v p is given by 



-Atot 
Btot 



(12) 



where 

A tot = A(0, 1, V&) - A(0, 1, Va) + A(as, r, <j> b ) - A(x, r,<j> a ) , 
Btot = B(0,l,ip b )-B(0,l,il> a )+B(x,r,cj> b )-B(x,r,ct> a ) 
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